(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     10986,        282]
NotebookOptionsPosition[     10453,        264]
NotebookOutlinePosition[     10808,        280]
CellTagsIndexPosition[     10765,        277]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"F", "[", 
   RowBox[{"f_", ",", "R2_", ",", "r_"}], "]"}], ":=", " ", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"\[Omega]", "=", 
       RowBox[{"2", "Pi", " ", "f"}]}], ",", 
      RowBox[{"RA", "=", "1"}], ",", 
      RowBox[{"R0", "=", "100"}], ",", 
      RowBox[{"R3", "=", 
       RowBox[{"10000", "/", "R2"}]}], ",", 
      RowBox[{"C", "=", 
       RowBox[{"10", "^", 
        RowBox[{"-", "7"}]}]}], ",", 
      RowBox[{"L", "=", "0.001"}], ",", 
      RowBox[{"U", "=", "1.78"}]}], "}"}], ",", 
    RowBox[{
     RowBox[{"NIntegrate", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Re", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"U", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "r"}], "-", "R2", "-", "RA", "-", 
                RowBox[{"\[ImaginaryI]", " ", "L", " ", "\[Omega]"}], "-", 
                RowBox[{
                "\[ImaginaryI]", " ", "C", " ", "r", " ", "R0", " ", 
                 "\[Omega]"}], "-", 
                RowBox[{
                "\[ImaginaryI]", " ", "C", " ", "R0", " ", "R2", " ", 
                 "\[Omega]"}], "-", 
                RowBox[{
                "\[ImaginaryI]", " ", "C", " ", "R0", " ", "RA", " ", 
                 "\[Omega]"}], "-", 
                RowBox[{
                "\[ImaginaryI]", " ", "C", " ", "R2", " ", "RA", " ", 
                 "\[Omega]"}], "+", 
                RowBox[{"C", " ", "L", " ", "R0", " ", 
                 SuperscriptBox["\[Omega]", "2"]}]}], ")"}]}], ")"}], "/", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "r"}], " ", "R2"}], "-", 
              RowBox[{"r", " ", "R3"}], "-", 
              RowBox[{"R2", " ", "R3"}], "-", 
              RowBox[{"r", " ", "RA"}], "-", 
              RowBox[{"R3", " ", "RA"}], "-", 
              RowBox[{"\[ImaginaryI]", " ", "L", " ", "R2", " ", "\[Omega]"}],
               "-", 
              RowBox[{
              "\[ImaginaryI]", " ", "C", " ", "r", " ", "R0", " ", "R2", " ", 
               "\[Omega]"}], "-", 
              RowBox[{"\[ImaginaryI]", " ", "L", " ", "R3", " ", "\[Omega]"}],
               "-", 
              RowBox[{
              "\[ImaginaryI]", " ", "C", " ", "r", " ", "R0", " ", "R3", " ", 
               "\[Omega]"}], "-", 
              RowBox[{
              "\[ImaginaryI]", " ", "C", " ", "r", " ", "R2", " ", "R3", " ", 
               "\[Omega]"}], "-", 
              RowBox[{
              "\[ImaginaryI]", " ", "C", " ", "R0", " ", "R2", " ", "R3", " ",
                "\[Omega]"}], "-", 
              RowBox[{"\[ImaginaryI]", " ", "L", " ", "RA", " ", "\[Omega]"}],
               "-", 
              RowBox[{
              "\[ImaginaryI]", " ", "C", " ", "r", " ", "R0", " ", "RA", " ", 
               "\[Omega]"}], "-", 
              RowBox[{
              "\[ImaginaryI]", " ", "C", " ", "r", " ", "R2", " ", "RA", " ", 
               "\[Omega]"}], "-", 
              RowBox[{
              "\[ImaginaryI]", " ", "C", " ", "R0", " ", "R3", " ", "RA", " ",
                "\[Omega]"}], "-", 
              RowBox[{
              "\[ImaginaryI]", " ", "C", " ", "R2", " ", "R3", " ", "RA", " ",
                "\[Omega]"}], "+", 
              RowBox[{"C", " ", "L", " ", "R0", " ", "R2", " ", 
               SuperscriptBox["\[Omega]", "2"]}], "+", 
              RowBox[{"C", " ", "L", " ", "R0", " ", "R3", " ", 
               SuperscriptBox["\[Omega]", "2"]}], "+", 
              RowBox[{"C", " ", "L", " ", "R2", " ", "R3", " ", 
               SuperscriptBox["\[Omega]", "2"]}], "+", 
              RowBox[{"C", " ", "L", " ", "R0", " ", "RA", " ", 
               SuperscriptBox["\[Omega]", "2"]}], "+", 
              RowBox[{"C", " ", "L", " ", "R2", " ", "RA", " ", 
               SuperscriptBox["\[Omega]", "2"]}]}], ")"}]}], "-", 
           FractionBox[
            RowBox[{"R2", " ", "U", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "R2"}], " ", 
                RowBox[{"(", 
                 RowBox[{"R0", "+", "R3", "+", "RA", "-", 
                  FractionBox["\[ImaginaryI]", 
                   RowBox[{"C", " ", "\[Omega]"}]]}], ")"}]}], "+", 
               RowBox[{"RA", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "R0"}], "+", 
                  FractionBox["\[ImaginaryI]", 
                   RowBox[{"C", " ", "\[Omega]"}]]}], ")"}]}]}], ")"}]}], 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"R2", " ", 
                 RowBox[{"(", 
                  RowBox[{"R0", "+", "R3", "+", "RA", "-", 
                   FractionBox["\[ImaginaryI]", 
                    RowBox[{"C", " ", "\[Omega]"}]]}], ")"}]}], "-", 
                RowBox[{"RA", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "R0"}], "+", 
                   FractionBox["\[ImaginaryI]", 
                    RowBox[{"C", " ", "\[Omega]"}]]}], ")"}]}]}], ")"}], " ", 
              RowBox[{"(", 
               RowBox[{
                SuperscriptBox["R2", "2"], "-", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"R0", "+", "R2", "-", 
                   FractionBox["\[ImaginaryI]", 
                    RowBox[{"C", " ", "\[Omega]"}]]}], ")"}], " ", 
                 RowBox[{"(", 
                  RowBox[{"r", "+", "R2", "+", "RA", "+", 
                   RowBox[{"\[ImaginaryI]", " ", "L", " ", "\[Omega]"}]}], 
                  ")"}]}]}], ")"}]}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"RA", " ", 
                 RowBox[{"(", 
                  RowBox[{"R0", "+", "R2", "-", 
                   FractionBox["\[ImaginaryI]", 
                    RowBox[{"C", " ", "\[Omega]"}]]}], ")"}]}], "-", 
                RowBox[{"R2", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "R0"}], "+", 
                   FractionBox["\[ImaginaryI]", 
                    RowBox[{"C", " ", "\[Omega]"}]]}], ")"}]}]}], ")"}], " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"R2", " ", "RA"}], "-", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "R0"}], "+", 
                   FractionBox["\[ImaginaryI]", 
                    RowBox[{"C", " ", "\[Omega]"}]]}], ")"}], " ", 
                 RowBox[{"(", 
                  RowBox[{"r", "+", "R2", "+", "RA", "+", 
                   RowBox[{"\[ImaginaryI]", " ", "L", " ", "\[Omega]"}]}], 
                  ")"}]}]}], ")"}]}]}]]}], "]"}], ")"}], "^", "2"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "0", ",", 
         RowBox[{"2", 
          RowBox[{"Pi", "/", "\[Omega]"}]}]}], "}"}]}], "]"}], "/", 
     RowBox[{"(", 
      RowBox[{"2", 
       RowBox[{"Pi", "/", "\[Omega]"}]}], ")"}]}]}], "]"}]}]], "Input",
 CellLabel->
  "In[922]:=",ExpressionUUID->"1179655b-cf84-493a-8b5d-914899bf95a7"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"f", ",", 
      RowBox[{"F", "[", 
       RowBox[{"f", ",", "100", ",", "0.01"}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"f", ",", "500", ",", "1000", ",", "10"}], "}"}]}], "]"}], 
  "]"}]], "Input",
 CellLabel->
  "In[924]:=",ExpressionUUID->"a39fafde-f32c-4840-842f-358ef0272384"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.011000000000000001`],
      AbsoluteThickness[1.6], PointBox[CompressedData["
1:eJwtz2tIk3EUx/FlNqxZzpQUzcvS0paat+yClx+RpSzJXKkhgYUamORmaBjB
ns25Z3elQl+U5DQVQ8EwjchEC+mi4AjKRKklVnQRb4TUNCv/58Dh8IHz5is5
X5pV6CIQCOL+7f+7NlDh6mrElOzGXNKaHSq0Bd2LjPBZYAaH4qHC+ZXuReZa
Dk8dx9L6i3/SP4eCmhrrpcQl5mg1JB/Scr2lv5g5NUqWr194FO9ktquhGhhs
unh2hTlYg+LYrE1JLavMCg26bXljNuG65DUPaKB0/nl9XObCLK5CwEwzF2Vb
z5xfhf3tWz6OuW1g7qpCo23ky4RGyCzQ4qh7T/mohxtzphYZbjny9K6NzI1a
1Enle6YKRMxzWhxRNmnjozYzoxr14S/7Jt09mGursSx1cn2cmNlRDbHQ0HpK
58kcrUOg10QnX7+VmdOhyCI63fbAi9mug9KzJEH43ps5mEfpdNHQktc2ZgWP
fXdzZPZsH+rn8XDmzPR0qy/16yE/KL5mcvWjfj0U/fmv/BT+1K9H8rA089vX
7dRvQM+dBXlWWSD1G1CReXj2jWsw9RuwyL/IC39HnjMgTBJ0P6NbQv1GjBvH
zZfrdlC/Ef63Rz4/V4dQvxHpz2QeleWh1G/Cod/nOnuUO6nfBOOk8ElHxS7q
NyGg94TFog2jfjMmdA39/K1w6jejsHcIw493U78Zda0HTuo/SanfgoS8kMpZ
3wjqt2BeaDXFZEdSvwXOlCvNqoYo6rdCEvJdIprdS/1WpJSJ3mpaoqnfisTl
jpv23BjqtyIutN3Be8ZSfw1Scxd/DI7GJv8FN7LzZQ==
      "]]}, {}}, {}, {}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{492.18750000000006`, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{500., 1000.}, {0, 4.379905951126329*^-9}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.7824741183615556`*^9, 3.7824744855542383`*^9},
 CellLabel->
  "Out[924]=",ExpressionUUID->"70795c31-bca7-4611-90ed-f4156c032bf3"]
}, Open  ]]
},
WindowSize->{1264, 637},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"11.3 for Microsoft Windows (64-bit) (2018\:5e743\:670828\
\:65e5)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 7287, 173, 409, "Input",ExpressionUUID->"1179655b-cf84-493a-8b5d-914899bf95a7"],
Cell[CellGroupData[{
Cell[7870, 197, 406, 12, 45, "Input",ExpressionUUID->"a39fafde-f32c-4840-842f-358ef0272384"],
Cell[8279, 211, 2158, 50, 223, "Output",ExpressionUUID->"70795c31-bca7-4611-90ed-f4156c032bf3"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)

